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Equivalent width measurements for rapid line variability in atomic spectral lines are degraded by increasing error bars 
with shorter exposure times. We derive an expression for the error of the line equivalent width o(W\) with respect to pure 
photon noise statistics and provide a correction value for previous calculations. 
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1 Introduction 

Error bars of line equivalent widths W\ have been the sub- 
ject of a number of investigations of high temporal and 
spectral resolution observations with exposures of the or- 
der of minutes or even seconds, leading to large errors. A 
first expression for the calculation of cr(W\) for pure pho- 
ton noise statistics was estimated by Williams et al. (1974) 
for their rapid line variability observations of Ap and Bp 
stars. This was then supplemented by Young (1974) for er- 
rors due to scintillation. Lacy (1977) published an improve- 
ment of Williams et al. and compared pure photon noise 
statistics with the case including scintillation for his Ha 
and H/3 rapid line investigations of massive stars. Guided 
by his work Chalabaev & Maillard (1983) again derived an 
expression for the error from pure photon noise statistics in 
an appendix of their paper on rapid spectral variability of 
Be stars. However, their derivation of a(W\) shows some 
inconsistencies which we attempt to resolve in this paper. 

2 The equivalent width of a spectroscopic 
line 

In practice, the measurement of noise within a line can be a 
difficult task, because photon noise is superposed with stel- 
lar variations. Hence, it is necessary to estimate an expres- 
sion for the equivalent width which separates these two pa- 
rameters and then to find an expression for its error a(W\). 
We start with the standard definition of the equivalent width: 
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with F C (X) the flux in the continuum, F(X) the flux in the 
line at the wavelength A and F(X) — F C (A) for A > A2 and 
A < Ai. In a first step we integrate eq. Q separately and we 
obtain 



W\ = AA 



F(X) 



dX. 



with the wavelength interval AA = 
the mean value theorem we obtain 



W\ = AA • 1 



F(X) 
F C (X) 



Ai. By applying 



(2) 



with ( pTxj ) the normalized average flux within AA. Equa- 



tion 10 can be written as 
It A ~ A A • \l- = 



(3) 



This can be shown by applying a discrete approximation 
of eq. l|3}- We use the arithmetic mean and substitute the 
line and continuum fluxes by their average values plus their 
deviations AF, and AF r .. 
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M represents the number of pixels within the line. Note that 
AFi within the line consists not only of noise but also of line 
information and can be quite large in contrast to AF Ci . 



F c M ^(l + ^SL) 



In case of sufficient S/N within the continuum we have 
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However, by definition 
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Hence, the equivalent width can be estimated from the 
wavelength interval and the average intensities in the line 
and the continuum. 

3 The error of the equivalent width 

The intensities F(\) and F C (X) are influenced by certain 
statistical errors. If we keep in mind that F c is generally 
estimated outside the line region and interpolated across the 
line where the line flux is measured and if we assume that 
the errors are not correlated we are able to estimate their 
standard deviations separately. Following the principal of 
error propagation we expand eq. in a Taylor series 
dW x , „ —. dW x 



W x = W{F, F c ) 



-=-(F - 

OF y 



F) 



oF r 



where F and F c are random variables. 

The variance Var(Wx) = a 2 (Wx) of the expansion is 
n 2 



<J 2 (W X ) 



1M> 



dF c 



~] 2 



(4) 



with (t(F) and <r(F c ) the standard deviations in the line and 
the continuum, respectively. By using eq. l|3} the two partial 
derivatives are 
dW x _AA 

with AA = M ■ hx and the spectral dispersion hx and 
W x ). 



^ = J=(AA 
dF n F, 



3.1 Weak lines 

For weak lines the depth of the line is negligible and 

F . T r 



a(F) 



S/N S/N 
with the signal-to-noise ratio S/N in the undisturbed contin- 
uum. With these terms we can write 
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(5) 



where the first term is the photometric uncertainty and the 
second term the uncertainty of the continiuum estimation 
over the line. According to eq. both terms are identical 
and we have 

(AA - Wx) 



a(Wx) =V2> 



S/N 



This result is easy to understand if we realize that the 
line strength for weak lines will be of the order of the noise 
itself. 



3.2 Low- and high-flux lines 

In the case of low (absorption) and high (emission) flux 
lines we have to use the corresponding poisson-statistic 




a(Wx) = 



(AA — Wx) 
S/N ' 



(7) 



This result seems to be quite simple. However, we 
should not expect a complex solution for the error bars of 
equivalent widths because of a relatively simple definition 
of W. In both cases we now can obtain the standard devia- 
tions with the measurable parameters S/N, AA, F and F c . 
In addition, eq. Q represents the general error of the line 
equivalent width and for the case F w F c we obtain the 
result for weak lines again. 



4 A comparison with former calculations 

Interestingly our result in eq. (0 differs from eq. (A9) of 
Chalabaev & Maillard (C&M) by a factor M, the number 
of pixels corresponding to the interval Ai < A < A2, in the 
first summand. The reason for this difference can be found 
in their expansion C&M (A3): 

- C&M (A3) is a multidimensional expansion for the 
equivalent width with average values Fj for all pixel 
fluxes Fj (the pixel fluxes Fj define random variables). 
However, it is not clear why the expansion over M 
fluxes Fj should be necessary. 

- Although there is no principal difference between the 
line and the continuum the continuum part is not con- 
sequently developed as it is done for the line part. They 
also say that "the value (S/N) characterizes the spread of 
read-outs on the continuum of an individual spectrum, 
while the quantity a(F c ) is the total uncertainty of the 
continuum level determination, and these two quantities 
may be related in a non-simple way". However, we as- 
sume that the uncertainty of the continuum is simply 
defined by the noise within the continuum and therefore 
the above distinction is not necessary. 

- Another difference between C&M and our approach can 
be found in their eq. (A7) which gives S/N — ^rjhj - 

However, the correct expression is S/N — (which 



is ,1 , in C&M notation). 



(6) 4.1 The VN-rult 



Generally, for the variance of a random variable x we have 
Var(x) = x 2 — x 2 with x the expectancy value of x. Re- 
spectively, we get Var(k ■ x) — k 2 ■ Var(x) with k = con- 
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stant. In particular, the mean value of x of N random vari- 
ables xi,...,xn1sx= jj ■ J2iLi x i and hence 



Var(x) = Var 



' 1 N \ 1 N 

8=1 / 1=1 



Note that x represents a random variable, as well. If all 
variances Var(xi) are identical the \/N-m\e gives cr(x) = 
a(x)/VN. 

In their approach C&M develop W\ in a multi- 
dimensional expansion 



9F C 



(Fc-i^c). 



We interpret the C&M notation of Fj and F c as follows: 

Fi = h- T.?=i Fi and T c = i • Eli F Ci with A> the 
number of pixel fluxes for the determination of the noise. 



Using 



dFj 



F c -M 



they get 



M 



F„ M 



dF c 



By calculating the variance of W\ one obtains 



Var(W x ) = f^=p-) " ■ Vor [ -1 • - Fj) 



aft 



y a r(F c -i^ c ). 



In their approach C&M now assume that both average 
values F and F c are constant (although this assumption is 
just an approximation for F by definition) and hence one 
gets 



yor(F c 



Now it is clear that C&M's goal was to estimate the vari- 
ance within the first summand with the variance over the 
averages F = -A= • Y^jLi Fj - The v^V-rule gives 
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and hence for the weak-line case 
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Fig. 1 Principle of the variance and average determination 
for a single spectrum and a number of spectra. The figure 
shows measured spectra with M pixels each. 



This is C&M (A9) and the first summand is smaller than 
our result in eq. (|5} by the factor M. For clarification we 
show the principal procedures in Figure^ Every measured 
single spectrum delivers a separate variance Var{F\). By 
obtaining N spectra, the variance for the wavelength-aver- 
aged spectrum is Var(F\) = jj-Var (F\). But considering 
the variance of all spectral averages, as is done by C&M, 
we get Var(F) — -h ■ Var{F\ ) independent of the number 
of measured spectra. Because of the missing sum within the 
continuum part of C&M (A3) the factor 1/M does not appear 
in both sums. 

If we want to estimate cr(W\) properly by means of 
C&M (A9), depending on the S/N over the whole line, we 
have to scale up the measured noise by the factor \AM. 
Therefore the corrected formula C&M (A7) is S/N = \f~M ■ 

F ■ 

a ip.\ ■ By using this correction Chalabaev & Maillard would 
have found the same result for a(W\) as we do. However, 
we note that this correction is enforced by the series expan- 
sion of C&M. Their procedure, however, has been used in a 
series of papers (more than 40 until now) using expression 
(A9) of Chalabaev & Maillard. 

Now it is easy to estimate the correction factor for the 
hitherto estimated error bars of W\ . For the normal case that 
both parts cr(F) and <j(F c ) are identical we find 



2M 
M + l 



C&M- 



For broad lines we have M + 1 w M and hence a {W\) w 
v2 • <j{W\)cSzM where the subindex C&M indicates the 
values obtained by Chalabaev & Maillard. In addition, if we 
assume that both errors (line and continuum) are different 
(as done by C&M) our result can differ even more. 

In fact, the above eq. Q for a 2 (W\) does not represent 
the exact variance of W but is just an approximation, which 
is especially important in the case of small basic popula- 
tions (i.e. the number of pixel fluxes for the determination 
of noise). The exact definition of the unbiased variance of 
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W\ is 



o 2 {Wx) 



N ■ (M — 1) f dW x 



N-M-l ' V dF 
M-(N-1) fdW 
N-M-l ' \dK 




(8) 



(Rreyszig 1967) with N the number of pixel values for the 
estimation of S/N and M the number of pixels within the 
line. This equation transforms into eq. @ if N and M be- 
come sufficiently large. If the number of pixels for the esti- 
mation of S/N and/or the line width are too small, the error 
of equivalent width will be overestimated by using eq. @. 
This is especially important for low -resolution observations. 
For example, if the investigated line covers only 10 pixel 
and S/N is estimated over about 50-100 pixels, the true er- 
ror would be overestimated by about 3%. 

4.2 Scrutinizing previous results 

Former conclusions regarding W\ have to be carefully scru- 
tinized due to increased error bars with respect to the result 
of C&M, especially if the interpretation of the data depends 
on results at the detection limit. For clarification, we give 
two prominent examples: 

- Searching for rapid line variability in the Be star BD 
+371 160, Rossi et al. (1991) tried to find a W\ correla- 
tion between Ha and H(3 (see their Fig. 4). Their con- 
clusion of "...spectral line variability on a time scale 
of 300 s for this Be star is unsustainable, especially for 



- Investigating periodic variabilities in UV lines of the B 
star £Cas, Neiner et al. (2003) claimed to see minima 
and maxima for W\ of the N V 1240 doublet (see their 
Fig. 2). This is not valid with respect to our correction 
factor of v2. 
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